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ABSTRACT 

The Planck and Herschel missions are currently measuring the far-infrared to millimeter emission of dust, which combined with 
existing IR data, will for the first time provide the full spectral energy distribution (SED) of the galactic interstellar medium dust 
emission, from the mid-IR to the mm range, with an unprecedented sensitivity and down to spatial scales ~ 30". Such a global SED 
will allow a systematic study of the dust evolution processes (e.g. grain growth or fragmentation) that directly affect the SED because 
they redistribute the dust mass among the observed grain sizes. The dust SED is also affected by variations of the radiation field 
intensity. Here we present a versatile numerical tool, DustEM , that predicts the emission and extinction of dust grains given their size 
£1^ distribution and their optical and thermal properties. In order to model dust evolution, DustEM has been designed to deal with a variety 

of grain types, structures and size distributions and to be able to easily include new dust physics. We use DustEM to model the dust 
SED and extinction in the diffuse interstellar medium at high-galactic latitude (DHGL), a natural reference SED that will allow us to 
study dust evolution. We present a coherent set of observations for the DHGL SED, which has been obtained by correlating the IR 
and HI-21 cm data. The dust components in our DHGL model are (i) polycyclic aromatic hydrocarbons, (ii) amorphous carbon and 
(iii) amorphous silicates. We use amorphous carbon dust, rather than graphite, because it better explains the observed high abundances 
of gas-phase carbon in shocked regions of the interstellar medium. Using the DustEM model, we illustrate how, in the optically thin 
limit, the IRAS/Planck HFI (and likewise Spitzer/Herschel for smaller spatial scales) photometric band ratios of the dust SED can 
disentangle the influence of the exciting radiation field intensity and constrain the abundance of small grains (a 5> 10 nm) relative to 
' the larger grains. We also discuss the contributions of the different grain populations to the IRAS, Planck (and similarly to Herschel) 

channels. Such information is required to enable a study of the evolution of dust as well as to systematically extract the dust thermal 
[ — , emission from CMB data and to analyze the emission in the Planck polarized channels. The DustEM code described in this paper is 

■ publically available. 

f"^ ■ Key words. Radiation mechanisms: thermal - Methods: numerical - ISM: dust, extinction - Infrared: ISM 

1 . Introduction sensitivity and angular resolution. Complementing these obser- 

vations with the existing IR data from ISO, Spitzer and IRAS, 
. ■ r-j , Dust plays a key role in the physics (e.g. heating of the gas, we win soon haye at our disposal the global dufjt SErx where 

^ ■ coupling to the magnetic field) and chemistry (formation of H 2 , the em i ss i n from grains of all sizes is observed. Such data will 

H ; shielding of molecules from dissociative radiation) of the inter- allow us t0 undertake system atic studies of the key dust evolu- 

, stellar medium (ISM). Heated by stellar photons, dust grains ra- tion processes in the ISM (e . g . grain growth and fragmentation), 

diate away the absorbed energy by emission in the near-IR to which are primarily reflected in the variation of the dust size dis- 

mm range. Dust emission can thus be used as a tracer of the tribution andi in particu i ari in the ratio of the abundance small 

radiation field intensity and, hence, of star formation activity. grains (a <; iq nm ) t0 i^ge grains ( e ,g, fctepnik et a . 2003; 

Assuming a constant dust abundance, the far-IR to mm dust p ara dis^taL| ||o09a; Ra pacioli et alJ I2005I: iBerne et all I2007I) . 



o 



emission is also used to derive the total column density along Dust eV olution encompasses many complex processes and its 

a line of sight and to provide mass estimates. The impact of dust impact may be studied using the dust SED for the diffuse inter . 

on the ISM and the use of its emission as a tracer of the local stellar me di U m at high galactic latitude (DHGL) as a reference, 

conditions depends on the dust properties and abundances. It is We present here a coheren t data set for the DHGL and the corre- 

therefore of major importance to understand dust properties and sponding dust mode i use d to match these observations, 
their evolution throughout the ISM. 

The instruments onboard the Herschel and Planck satel- A suitable dust model has to, ideally, explain all of the 

lites are currently observing the dust spectral energy distribu- available observational constraints, these are primarily: (i) the 

tion (SED) from the FIR to the millimetre with unprecedented extinction curve, (ii) the emission brightness spectrum vl v or 
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Table 1. Dust SED, vl v , and its uncertainty, cr{vl v ), in unit of 10~ 30 W sr~ l H l measured for the Diffuse High Galactic Latitude 
(DHGL) medium. The DIRBE and WMAP values were obtained by correlation of the measured intensities with H I emission data 
and corrected to account for the contributions of the diffuse ionized medium and the diffuse molecular gas (see § 14. 11 1. The Herschel 
PACS, SPIRE and Planck/HFI values were derived by integrating the FIRAS spectrum (see § 14. Il l considering the instrument trans- 
mission and color correction. 
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SED (Spectral Energy Distribution), (iii) the elemental deple- 
tions, (iv) the spectral dependence and efficiency of polarization 
in emission and extinction and finally, (v) the plausibility of a 
given dust component in view of its likely formation and de- 
struction in the ISM. Early models focused on explaining the 
extinction (some of them also describing polarization by prefer- 
entia l extinction) and used silicates and carbonaceous material 



e.g. iMafhis eTaLlll977t iDraine & Led 119841; iKim et alj [19941 



Kim & Martinlll995t ILi & Greenberg|ll997l) . In the 1980's, the 



spectroscopy of dust emission in the IR lead to the introduc- 
tion of a population of aromatic hydrocarbons called polycyclic 
aromatic hydrocarbons (PAHs). However, current models aim 
to consistently rep roduce the observe d extinction and emission 
(IDesert et al.l 19901 hereafter DBP90) (ISiebenmorgen & Kruegell 
1992tlDwek etalJ|l997t ILi & DraineH2001at IZubko et al.l l2004: 
Draine & Lill2007l hereafter DL07 



as well as the associated po- 
Further details of the dust 



.'7) 

larization dDraine & Fraissell2009l) . 
modelling process can be found in rec e nt reviews on interstella r 
dust, such as those bv lDraind d2003al) ; ILi & Greenbergl d2003l) : 
iDraind (120041) . We emphasize here that all of the different dust 
models, developed to date, have been de signed using essentially 
the same observa t ional constraints ( e.g. ILi & Greenbergl [ 1997; 
IZubko et al.ll2004lDraine & Lill2007h . 

Carbonaceous dust is often considered to be in the form 
of crystalline graphite because it can explain the 217 nm bump 
in the observed extinction curve. The question of the origin 
of this feature is however still open because graphite can- 
not explain the band profile variability o bserved in the ISM 
(Fitzrjatrick & Massa 2007; Draine 2003a). In fact, throughout 
the dust lifecycle, from dense molecular cores to the diffuse ISM, 
the grain populations are the result of complex, non-equilibrium 
evolutionary processes and it appears natural to consider that 
carbon dust is amorphous, as observations have clearly demon- 
strated for silicate gr ains (e.g. lLi & Drainel2001bl;lKemper et alj 
l2005HLi el al. 12007). We follow this train of thought in this work 
and examine the possibility that all carbonaceous grains (except 
PAH) are amorphous and partially hydrogenated in the interstel- 
lar medium. 

The paper is organized as follows: we first present our model 
DustEM and the developed dust SED-fitting tool (Sec©. We 
then discuss why amorphous carbon is the most likely car- 
bonaceous dust grain component in the ISM (Sect[3]l. Section|4] 
presents a coherent data set for the DHGL SED and the asso- 
ciated dust model that serves as a reference SED for the eval- 
uation of the nature and efficiency of dust evolution processes. 
In Sect.[5] we show how changes the dust abundances and sizes 
can be traced from their global SED, thus allowing a study of 



dust evolution. The conclusions and perspectives of this work 
are drawn together in Sect|6] 



2. The DustEM model capabilities 

Accounting for the emission from dust, for which the proper- 
ties evolve in the ISM, requires a model that easily allows the 
dust properties to be changed and mixed, and where new dust 
physics can easily be incorporated and tested. In this section we 
present a numerical tool, DustEM , which is designed to this end 
and that provides the extinction and emission for interstellar dust 
populations. Note that DustEM computes the local dust emissiv- 
ity, and as such, does not include any radiative transfer calcula- 
tion, so that it can be included as a source and extinction term in 
more complex applications explicitly dealing with energy trans- 
port through optically thick media. 

2.1. The DustEM model 

The DustEM model stems from that of DBP90 but the code has 
been completely re-written in fortran 95 in order t o achieve the 
goals described above. We use the formalism of IDesert et al.l 
(1986) to derive the grain temperature distribution dP/dT and 
DustEM then computes the dust SED and associated extinction 
for given dust types and size distributions. To correctly describe 
the dust emission at long wavelengths the original algorithm has 
been adapted to better cover the low temperature dP/dT region. 
Usin g an adaptative temp erature grid, we iteratively solve Eq. 25 
from Des ert et al.l ([1986) in the approximation where the grain 
cooling is fully continuoufl This topic will be discussed in de- 
tail in a forthcoming paper, along with its consequences for the 
emission of small grains and gas-grain interactions. Note that un- 
less explicitly specified in the input file, the dP/dT calculation is 
performed for all grain populations and sizes including those for 
which the thermal equilibrium approximation would apply. The 
principal input parameters for the dust model (e.g. grain types, 
dust-to-gas mass ratios) used to compute the required quantities 
are read from a control input file along with keywords (e.g. that 
specify the size distributions or the outputs). DustEM can handle 
an arbitrary number of grain types. The built-in size distributions 
dn/da (the number of grains of radius between a and a + da) are 
either power law or log-normal but any size distributions can be 
given in input data files. For each grain of specific mass den- 
sity p, the radius a is defined as that of an equivalent sphere of 



The grain cooling is taken to be as in Eq. 8 of IDesert et all Jl986h 
where the heating term is neglected and the cut-off energy is that of the 
hardest photons in the radiation field. 
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Table 2. DHGL dust model abundances and size distribution pa- 
rameters (see 34.2b . Y is the mass abundance per hydrogen for 
each dust component. fu m is the dust component mass as a frac- 
tion of the total dust mass. 
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mass m = p4na 3 /3. For a given dust type, the optical proper- 
ties are provided as a function of the grain radius a in the data 
files, namely: Q-files for the absorption and scattering efficien- 



cies, Q(a, A), and C-files for the heat capacity per unit volume, 
C(a, T The Q and C values are provided over the broadest 
possible size range (usually ~0.3 to 10 4 nm). Any grain type can 
be considered once these Q and C-files have been generated and 
are in the DustEM data file set. The requested dust sizes have to 
fit within those used in the above files (i.e. extrapolation is not 
allowed). 

A size-dependent weight factor (or mixing) read from a data 
file can be applied to any grain type while computing its ex- 
tinction and emission. This allows us to take into account a 
number of effects, e.g., the ionization of PAHs or more gen- 
erally the evolution of one grain type to another as a function 
of size. The format of the DustEM tabulated size distributions 
(and their elemental composition as a function of size, in case 
of composite grains) data files is shared with DustEV, the dust 
evolution model described in Guillet et al. (2010, in prepara- 
tion) and which follows the evolution in composition and size 
of charged dust grains as a result of coagulation and fragmenta- 
tion processes. The default outputs from DustEM are the emis- 
sivity per hydrogen atom (A^h = Nui + 2 A^h, ) for each grain 
type Anvlv/Nu in ergs -1 H 1 and the extinction cross-section 
per gram for each grain type (in cm 2 g _1 ). When required (and 
via keywords in the input control file) the grain temperature dis- 
tribution and the dust emissivity can be provided for each grain 
size. When required DustEM also allows the use of temperature- 
dependent dust emissivity (see AppendixIBl. 

The current development plan for DustEM consists of the in- 
clusion of (i) the spinning dust emission, (ii) the polarized ex- 
tinction/emission, and (iii) the ph ysics of amorphous solids at 
low temperatures ( Men vetafl 120071) . As for the mixing or the 
/^-correction, each process is activated with a keyword given in 
the input control file and the associated parameters or tabulated 
values are then given in a specific file. DustEM can also be used 
as a subroutine in radiative transfer models, such as the Meudon 
PDR codfl which solves for the gas state and dust emissivity. 
DustEM can also be use d as a subroutine of the 3D-continuum 
transfer model (CRTfl of Uuvela & Padoanl (120031) . 



2 Size-dependent heat capacities are required to describe PAH (C-to- 
H ratio) and composite grains. 

3 Available at http : //par . obspm . fr/PDRcode . html 

4 Available at 

http: //wiki .helsinki . £i/display/~mjuvela@helsinki . fi/CRT 
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Fig. 1. Mass distribution for the four dust types used in the 
DHGL model. 



2.2. The IDL wrapper 

The DustEM code is coupled to a wrapper written in IDL 
(Interactive Data Language). The wrapper allows the user to run 
DustEM iteratively from a code development and debugging en- 
vironment. In particular, it is possible through this interface, to 
provide a given SED, and the associated erro r bars, to be fitted 
by DustEM . The SED fitting uses the mpfit (iMarkwardfe OO ) 
IDL minimization routing, which is based on the Levenberg- 
Marquardt minimization method. This allows the user to itera- 
tively modify any of the model's free parameters, in order for 
the model prediction to converge on the provided SED. In this 
process the wrapper uses the DustEM computed spectrum to pre- 
dict the SED values that would be observed by astronomical in- 
struments, iteratively taking into account the appropriate color 
correction for wide filters and the flux conventions used by the 
instruments. The currently recognized settings include any fil- 
ter used by the IRAS, COBE (DIRBE, DMR), Spitzer (IRAC, 
MIPS), Planck (HFI, LFI) and Herschel (PACS, SPIRE) satel- 
lite instruments, as well as that of the Archeops and Pronaos 
balloon-borne experiments. The target SED can also include ar- 
bitrary spectroscopic data points, (such as for the COBE-FIRAS 
or Spitzer-IRS observations) for which no color correction is ap- 
plied. The fitting of an SED to a reasonable accuracy through 
this interface generally requires a few minutes of a single CPU. 
The resulting best fit parameters are provided with their 1-sigma 
error bars, as derived from the mpfit routine. The wrapper also 
uses the capability of mpfit to tie various free parameters through 
a user specified function. It also includes the concept of user- 
written subroutines (e.g. a subroutine to manage the PAH ion- 
ization fraction). The wrapper can also be used to generate 
pre-calculated SED tables for the instruments mentioned above, 
which can then be used to perform regression for a limited set 
of free parameters. This method is well adapted to problems 



To be found at http : //purl . com/net/mpf it 
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10 100 1000 

\ (jum) 

Fig. 2. Dust emission for the DHGL medium. Grey symbols and curves indicate the observed emission spectrum (see S I4.1I ) for 
N H = 10 20 H cnT 2 . The mid-IR (-5-15 pm) and far-IR (~ 100 - 1000 /mi) spectra are from ISOCAM/CVF (ISO satellite) and 
FIRAS (COBE satellite), respectively. The triangle at 3.3 jt/m is a narrow band measurment from AROME balloon experiment. 
Squares are the photometric measurments from DIRBE (COBE). Black lines are the model output and black squares the modeled 
DIRBE points taking into account instrumental transmission and color corrections. 



where the direct fitting of individual SEDs over many image pix- 
els would require too much CPU time. 

Past applications of the DustEM wrapper to fit SEDs have 
included deriving the InterStellar Radiation Field (ISRF) scal- 
ing factor and the dust mass abundances (Y) for the various dust 
grain types, as well as constraining the size distribution of very 
small grains (VSGs), based o n SEDs including a combination 
of IRA S and Spitzer data (e.g. lBernard et al . 2008; Pa radis et alj 
12009b!) and Herschel data (Bernard et al. 2010). Future devel- 
opments of the wrapper will include the possibility to simulta- 
neously fit extinction and emission and to follow the planned 
evolution of the fortran code. 

DustEM is developed within the framework of a collab- 
oration between IAS and CESR. The fortran code, the IDL 
wrapper and the associated data files can be downloaded from 
http : //www . ias . u-psud . fr/DUSTEM. 



3. The nature of interstellar carbon grains 

We discuss here the plausible nature of carbon dust grains 
in the ISM. Carbonaceous dust is formed in the ejecta of 
evolved stars but, seemingly, must also form directl y in the ISM 
(IZhukovska et al.ll2008t |Prainell2009i3: lJonesll2009h . The small- 
est carbonaceous species contain a few tens of C-atoms and are 
assigned to PAH-like species due to their spectral similarity to 



the observed mid-IR emission bands. The nature of the larger 
carbonaceous grains is less clear, in particular, because of the 
lack of specific spectroscopic signatures. Graphite was long-ago 
proposed as an interstellar dust component because its optical 
anisotropy could nat urally explain the polarizat ion of starlight 
dimmed by dust (e.g. lCavrel & Schatzmanl fl954) and especially 
since it features a strong absorption band (due to n — » n* tran- 
sitions in sp 2 carbon) similar to the 217 nm bump observed in 
interstellar extinction dStecher & Donn|[l965l) . However, such a 
feature exists in many carbonaceous materials containing aro- 
matic rings and is not specific to graphit e. In addition, de- 
tailed modelling of the band profile by iDraine & Malhotral 
(11993b has shown that graphite grains cannot explain th e 
observed bandwid t h var iability (Fitzpatr ick & Massal 120071) . 
IDraine & Malhotral d 1993b suggest that this variability may be 
due to the presence of hydrogen or defects within graphite. 

In fact, observational evidence strongly suggests that carbon 
dust consists of hydrogenated amorphous carbon (HAC or a- 
C:H) mate rials that form a round evolved stars, metamorphose 



there (e.g., 



iJones et a! 



Goto et al. 2003) and further evolve in the ISM (e.g., 
199(1 iPino etal1l2008l: Ijonesll2009h . HAC material 



could also form directly into the ISM from the UV processing 
of ice mantles after lea ving the dense medium as suggested by 
iGreenberg et al J (119951) or via direct accretion of carbon and hy- 
drogen (e.g. lJones et al.ll 1990h . In quiescent regions of the dif- 
fuse ISM ~ 40% of the available carbon is observed to be in the 



4 



M. Compiegne et al.: The global dust SED: Tracing the nature and evolution of dust with DustEM 




100 



1000 
X (yitm) 



Fig. 3. Long wavelength emission of dust in the DHGL medium. 
Grey symbols and curves indicate the observed emission spec- 
trum (see § 14. Il l for Nh - 10 20 H ctrT 2 . The spectrum is a 
FIRAS (COBE satellite) measurment. Squares and diamonds 
are the photometric measurments from DIRBE (COBE) and 
WMAP, respectively. Black lines are the model output and black 
squares the modeled DIRBE points taking into account instru- 
mental transmission and color corrections. 



gas (e.g JCardelli et al1 [T996). In contrast, in regions of the ISM 
shocked to velocities of the order of 100 - 150 km s _1 as much 
as 80 — 100% of the carb o n atoms are observed to be in the gas 
phase dFrisch et alj |1999t IWeltv etal.ll200H iPodioet al.ll2006t 
ISlavinll2008l) . indicating that a large fraction of the carbon dust 
has been de stroyed. Graphite/amorphous c arbon dust p rocess- 
ing models dTielens et al.lll994H Jones et alJI 1994 1 19961) predict 
that, at most, only 15% of graphitic carbon dust should be eroded 
in such shocks (for silicate dust the observations and model re- 
sults agree rather well). This discrepancy for carbon dust im- 
plies that we ne ed to re-consider the nature of carbonaceous 
dust in the ISM. ISerra Dfaz-Cano & Jonesl d2008l) re-evaluated 
the processing of carbon grains in shocks, using the physical pa- 
rameters typical of a-C:H/HAC, rather than of graphite. Their 
newly-calculated sputtering yields for a-C:H/HAC indicate that 
it is significantly more susceptible to sputtering than graphite 
and they showed that 30 - 100% of these type of carbona- 
ceous grains are destroyed (= 60 - 100% carbon in the gas) 
in 100 - 200 km s _1 supernova-generated shock waves in the 
warm inter-cloud medium. These results, along with similar re- 
sults obtained for PAH evolution in shocks and in a hot gas 
dMicelotta et al . 20 10b, a), are in good agreement with the obser- 
vations of shoc ked regions of the ISM (e.g., 80-100% of the car- 
bon in the gas:lFrisch et alj |1999t IWeltv et alj|2002t iPodio et alj 



2006: ISlavinll2008l) 

Thus, even if all the carbonaceous dust formed around 
evolved stars were to be in graphitic form, an unlikely sce- 
nario^! would subsequently be modified, by ion irradi ation (e.g., 
Banhartl [l999l: iMennella et ai]|2003t lMennellal l2006) in shocks 
and by cosmic rays in the ISM, to an amorphous form im- 
planted with heteroatoms, principally H atoms. However, given 
the results for the processing of a-C:H/HAC dust and PAHs 
dSerra Dfaz-Cano & Jonesl l2008t iMicelotta et all l2010bllah . the 
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Fig. 4. The extinction curve (top) and albedo (bottom) for dust 
in t he DHGL medium . In the top panel, the grey diamonds show 
the iFitzpatrickl dl999l) extinction law for Ry = 3.1. Error bars 
come from the dispersion of the observed extinction curve. The 
bottom panel overlays the model albedos from DBP90 and DL07 
(grey lines) and from the DustEM model (black line), compared 
to the observational data points. 



ion and electron irradiation of carbon grains in shocks appears 
to be very destructive. Therefore, the carbon dust in the ISM 
is probably completely re-processed on short time-scales, rather 
than simply being modified by ion irradiation effects in shocks. 

The physical and chemical properties of a-C:H/HAC solids 
are sensitive to the external conditions (e.g., photon, ion or 
electron irradiation) and can undergo a process of 'aromatiza- 
tion' when heated or exposed to UV or ion irradiation (e.; 
iDulev et al.Hl 9891: lJoneslfl990l iJones etaDll99Ct iJonesI 12009 
Labo ratory experiments on interstellar carbon dust analogues 
(e.g jDartois et alj |2004) show that hydrogen-rich (> 50 atomic 
% H) HAC solids can explain the interstellar absorption bands 
at 3.4, 6.85 and 7.25 fim and that the thermal annealing of the 
material is accompanied by an increase in the aromatic carbon 
content (i.e., 'aromatization' ). The HA C model for hydrocar- 
bon grains in the ISM (e.g., |Jonesll2009l and references therein) 
would therefore predict that the larger grains, with tempera- 
tures in equilibrium with the local radiation field, should be 
rather hydrogen-rich because they do not undergo significant 
heating. In contrast, the smaller, stochastically-heated grains, 
will be converted to hydrogen-poorer, lower density, smaller 
bandgap, aromatic -rich materials. The photo-fragmentation of 
the smaller, aromatic -rich grains could then be an important 
source of molecular aromatic species such as PAH s and molec- 
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Fig. 5. Infrared extinction of dust in the DHGL medium. The 
grey circles show the observed extinction (Mathisl [l990h . Black 
lines are the model. 



ular fragments (e.g..lJoneslll990tlpulevl l2000t iPetrie et alj|2003t 
Petv et al.|[2005t I Jonesll2005L l2009l) . 

In addition to the PAHs, we assume that the interstellar 
carbonaceous grains are in the form of hydrogenated amor- 
phous carbons, collectively known as a-C:H or HAC, rather 
than graphite, and adopt the se materials as an analogue for 
interstellar carbon dust (e.g., Honesl 19901; Dartois etal. 2004; 



Dartois & Munoz-Carol 120071; iDartois et alj 120071; iPinoetall 
2008; lJonesll2009l) . For the optical properties of these amor- 
phous carbons, we here use the BE sample complex refractive in- 
dex da ta, derived from laboratory measurement by Zubk o et al.l 
(1996j), to be representative of sp 2 and sp 3 containing (hydro- 
genated) amorphous carbons (see Appendix lA.21 >. However, we 
point out that these materials can have a wide range of properties 
and these can evolve as a function of the local conditions. 



4. A dust model for the diffuse high galactic latitude 
medium 

We follow earlier studies (e.g. DPB90) in using observations 
of the diffuse interstellar medium in the Solar Neighborhood to 
constrain our dust model. The model fitted to these observations 
defines a reference SED with which to compare and characterize 
dust evolution as a function of the local interstellar environment 
in the Galaxy and in external galaxies. In this section, we de- 
scribe the observations used and the model that reproduces these 
observations. 

4.1. The observables 

A first observational constraint on the dust model comes from 
measurements of the dust extinction and scattering properties. 
For the ext i nction curve we use the data , for Ry = 3.1 , from 
iFitzpatrickl dl999l) in the UV-visible and iMathisi d 19901) in the 
infrared. We normalize the observations using the empirical de- 
termination of N H /E(B - V) = 5.8 x 10 21 H cm" 2 , obtained 
by correlating the reddening of stars with the HI and H2 column 
densities derived fro m Copernicus and FUSE far-UV spectro - 
scopic observations dBohlin et al.lll978"l iRachford et al.l 12002). 
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Fig. 6. The dust emission for our DHGL model compared to 
those of DL07(U=1,4p AH = 4.6 %) and DBP90 for N H = 10 20 
H cm 2 . For DBP90, the measured SED used as a reference 
was higher explaining the deviation at wavelengths larger than 
~ 20 yum. 



The dus t albed o is taken from the wor k of Lillie & Wittl dl976l) : 
lMorganldl980l) ; IChlewicki & Laureiisldl988l) . 

The dust model is also constrained by observations of the 
dust emission from near-IR to mm wavelengths. We use data at 
high Galactic latitudes from the COBE and WMAP space mis- 
sions todeterminejhe SED ofjhedtistemission. Several authors 
(e.g. Boulange r et al.l fl996t lArendt et aD [l998) have analyzed 
the COBE data, but, to our knowledge, no publication presents a 
coherent set of measurements obtained over the same sky area. 
To reduce the uncertainties associated with variations of the dust 
emission per Hydrogen atom across the high Galactic latitude 
sky, we have made our own determinations. We used the DIRBE 
and FIRAS data sets, corrected for zodiacal light, available at 
lambda . gsf c . nas a . gov. For WMAP , we used the seven-year 
temperature maps dJarosik et alJ feOlO). and subtracted the cos- 
mic microwave background (CMB) usin g the map obtained with 
the internal linear combination method dGold et all |20 1 0) . The 
free-free contribution to t he emission was subtracted using the 
model of iDickinson et al.l d2003l) . 

At wavelengths > 60 /itn up to the WMAP centimetric data 
we determined the interstellar emission per Hydrogen atom by 
correlating the dust emission with the HI line emission, Ihi, 
taken from the Leiden- Argentina-Bonn survey (Kalberla et al. 
l200l . We obtained a coherent set of measurements by using 
the same sky area, defined by the galactic latitude and HI emis- 

sion constraints \b\ > 15 and Ihi < 300 K km s , for all 
wavelengths. We discarded bright point sources, nearby molec- 
ular clouds in the Gould Belt, and the Magellanic clouds. The 
error-bars on each measurement were obtained by dividing the 
sky area into four sub-areas, and comparing results between sub- 
areas. These error-bars are dominated by systematic variations of 
the SED from one sub-area to the others. 

It is not possible measure the interstellar emission at near- 
and mid-IR wavelengths over the same sky area that we used 
for the longer wavelengths. At mid-IR wavelengths, and outside 
clean areas near the ecliptic poles, uncertainties in the zodiacal 
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Table 3. The elemental abundances (see § 14. ll for details). 





c 





Mg 


Si 


Fe 


Measurements 


°[X/10 b H] o 
"[X/10 6 H] F , G * 


269±33 


490±60 


40±4 


32±2 


32±3 


358±82 


445±156 


42.7±17.2 


39.9+13.1 


27.9±7.7 


c [X/10 b H] gas 


75±25 


319+14 


~0 


~0 


~0 


[X/10 b H] dust/O 


194±41 


171±62 


40+4 


32+2 


32±3 


[X/10 H]dust/F,G* 


283±86 


126±157 


42.7±17.2 


39.9+13.1 


27.9±7.7 


DHGL dust model 


[X/10 b H] PAH 


65 










[X/10 6 H] S amC 


14 










[X/10 6 H] L!lraC 


121 










[X/10 6 H] aSl , 




180 


45 


45 


45 


[X/10 b H] TO TAL 


200 


180 


45 


45 


45 



Notes. M fromlAsplund et alj f2009) 

{b) from lSofia & Mevd d20oTh 

(c) C estimate is from lDwek et all (1997) while the O estimate is from lMever etaD dl998h 



light removal are larger than the high galactic latitude interstellar 
emission. At near-IR wavelengths the stellar emission is domi- 
nant and its subtraction is an additional difficulty in measuring 
the interstellar emission. At wavelengths < 25 fim, we have thus 
used the colors, I v (A)/I v (10Q p.m), measured over distinct sky 
areas by lArendt et al. d!998l) (see their Table 4). We point out 
that their 140/1 00 fim and 240/1 00 fim colors, measured over a 
smaller sky area near the Galactic poles |b| > 45°, are 10 - 15% 
lower than our values. This difference reflects small variations of 
the far-IR colors across the high Galactic latitude sky. These spa- 
tial variations make the near and mid-IR parts of the SED more 
uncertain than the far-IR part, which we determined consistently 
over a single sky area. We have complemented the DIRBE data 
points with the narrow band (AA = 0.17 fim) photometric point 
at 3.3 /im from the AROME balloon observ atory measured to- 
wards the molecular ring dGiard et alJ ll994). The mid-IR spec- 
trum is the ISOCAM/CVF spectrum for a diffuse region at g alac- 
tic coordinates (28.6, +0.8) presented in Fla gev etall d2006l) and 
scaled to match the 12 fim emission per hydrogen. 

For optically thin emission the HI cut-off used in the correla- 
tion analysis corresponds to a column density of neutral atomic 
hydrogen of 5.5 x 10 20 ctrT 2 . Over this sky area, the interstellar 
emission arises primarily, but not exclusively, from dust within 
neutral atomic interstellar gas. Emission from dust associated 
with diffuse H2 and HII gas also need s to be taken into ac - 
count. Using the FUSE data reported by Gillmon et al. (2006), 
we estimated the fraction of H2 gas over the selected sky area 
to be ~ 3%. Io nized gas accounts for a larger fraction (~ 20% 
lRevnoldsl[T989]) of the hydrogen in the diffuse ISM. The U a 
emission from the diffuse HII gas has been shown to b e spatially 
correlated with HI gas emission (Rey nolds et alJ[l995l) . We as- 
sume that the dust emission per hydrogen atom is the same in 
the neutral and ionized components of the diffuse interstellar 
medium, and thus we scale the SED derived from the dust-HI 
correlation by 0.77, to obtain a dust emission per hydrogen atom 
that takes into account the contributions of the diffuse ionized 
medium and the diffuse molecular gas. The corrected SED val- 
ues and error-bars are listed in Table[TJ where the FIRAS spec- 
trum has been convolved with the Herschel and Planck transmis- 
sion for the relevant spectral bands. The full SED is plotted in 
Fig.|2]and Fig. [3] From 100 fim to 3 mm, the SED is well fit by a 
black body spectrum with a temperature of 17.5 K, and an emis- 



sivity oc v 2 for A < 450 yt/m and oc y 1 - 65 at longer wavelengths. 
The diffuse HII must also contribute to the extinction along dif- 
fuse lines of sights as suggested by the difference between in- 
tercloud an d cloud for the valu e of (N(HI) + N(H 2 ))/E(B - V) 
reported by Bohlin et alJ (1 19781) . Here we use the standard value 
of 5.8 X 10 21 ctrT 2 mag 1 also found bv lRachfordetalJ d2002l) 
towards "translucent lines of sight" (/# 2 ~ 0.5, 1 S= Ay ^ 3). 

Table[3]lists the measured elemental abundances and inferred 
dust elemental abundan ce relevant for our m odel. The solar 
abundances are those of lAsplundetaT] (l2009h while the F, G 
stars (young < 2Gyr F, G disk stars) abundances are those of 
[Sofia & Meyeij d2001h . We do not consider B stars that have 
lower metallicity, because of the possible segregation of refrac- 
tory elements during the star formati on due to se dimentation 
and/or ambipolar diffusion processes dSnowl 12000). Note that 
this could also be the case for lower mass stars (e.g. the sun, 
Asplund 2008), raising the question of the existence of a stellar 
standard that represents the overall ISM abu ndances. Th e inter- 
stellar gas phase abundan ces used are from Dwe k et alJ d 19971) 
for the carbon and from iMever et all d!998j) for the oxygen. 
The iDwek et al.l d 19971) value of 75+25 ppm for the gas phase 
C abundance is in goo d agreement with recent estimates by 
ISofia & Parvathil d20 09). We emphasis that, as argued bv lDraind 
d2009bl) . in view of the assumption that are made for dust prop- 
erties (especially density/porosity) and the uncertainties on the 
abundances, a model that departs from the measured elemental 
solid phase abundances by tens of percent should still be consid- 
ered as viable. 



4.2. The model 

To model the DHGL data set we use three dust components: 
(i) PAHs, (ii) hydrogenated amorphous carbon (hereafter amC) 
and (iii) amorphous silicates (hereafter aSil). The properties of 
each type are described in the AppendixlAl In the course of dust 
evolution, composite or core-mantle (silicate-carbon) grains are 
likely to appear that would lead to a similar polarization of both 
the 3.4 fim absorption ba nd of amorphous carbon and the 9.1 fim 
band of the silicates (see Adams on et al ] ll999HLi&Greenbergl 
2002). We do not consider such composite grains here because 
recent spectropolarimetric measurements by IChiar et al.l ([2006) 
indicate that the 3 A fim band is not polarized while the 9.1 fim 
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Fig. 7. The modeled dust emission per H atom for scaling factors for the the MMP83 ISRF of {/=0.1 to 10 4 . All spectra have 
been divided by U to emphasize the changes in the shapes of the SEDs, which are presented as I v to better indicate the photo- 
metric data points. The dust properties are those for DHGL model (see § I4.2| |. The left panel shows the Spitzer (triangles) and 
Herschel (squares) photometric points. The right panel shows the IRAS (triangles) and Planck/HFI (squares) photometric points. To 
emphasize excitation effects we have assumed a constant ionization fraction for the PAHs. 



band along the same sightline is polarized. Although such ob- 
servations should be pursued, these results may indicate that the 
carriers of the 3.4 and 9.7 fim features belong to physically sep- 
arated dust populations where the former are poorly elongated 
or not well aligned with the magnetic field while the latter are. 
Further constraints on this issue will be obtained from Planck 
measurements of polarized dust emission in the submm range. 

The size distribution and abundance of the dust components 
have been adjusted in order to reproduce both the extinction and 
emission. The required size distributions are shown in the form 
of their mass distributions in FigQ] The population of amor- 
phous carbon dust has been divided into small (SamC) and large 
(LamC) grains, but the overall size distribution of carbonaceous 
grains ( PAH and amC) is continuous . Following previous stud- 
ies (e.g. Wein gartner & Draind2001al) . the smallest grains (PAH 
and SamC) size distributions are assumed to have log-normal 
distributions (with Oq the centre radius and cr the width of the dis- 
tribution). The LamC and aSil size distributions follow a power 
law (oc a a ) starting at a m \ n and with an exponential cut-off of the 
fo rm e~K a " fl ')/ fl ']' - for a > a t (1 otherwise) at large sizes. As noted 
by iKim & Martini d 19951) . such a form in the cut-off is required 
to explain the observed polarization in the near-IR: the polariza- 
tion capabilities of the present dust model will be discussed in 
a future paper. The abundances and parameters for the size dis- 
tributions are summarized in Table[2] PAHs and SamC are small 
grains of radius less than ~ 10 nm while LamC and aSil cover a 
broad size range from a few nm to a few 100 nm. The mean size 
of PAHs (flo = 6.4 A) corresponds to 120 C atoms. Note that, 
contrary to carbonaceous dust, for which the PAH bands give 
a constraint on the abundance and size of very small particles, 
only an upper limit can be set regarding the abundance of very 



small silicate dust (see lLi & Draind2001bl) . In our model, we do 
not include very small silicate particles. 

Our model results are compared to the observational data in 
Figs. [2] [3] |4] and [5] The dust temperatures were c omputed us- 
ing th e InterStellar Radiation Field (ISRF) from Mathi s et al.l 
(119831) . hereafter MMP83, for D c = lOkpc, labelled "solar 
neighborhood". Between 5 and 16jum, we convolved the model 
SED with a Gaussian profile {A/AA ~ 40) to account for the 
ISOCAM/CVF spectral resolution. The LamC and aSil grains, 
mostly in thermal equilibrium with the ISRF, emit in the far-IR 
while PAHs and SamC grains undergo temperature fluctuations 
and thus contribute to the near- to mid-IR emission. The LamC 
grains dominate the observed flux at 100 and 140 jum since they 
are hotter than aSil thanks to their higher absorption efficiency 
at A ~ 1 fan (see Fig.lAlll where the MMP83 ISRF peaks. We 
note that the SamC grains make a significant contribution at 25 
and 60 fim. At A <: 250 fim, the aSil and LamC contribute almost 
equally to the SED down to the WMAP measurements that are 
well reproduced (see FigO denoting no significant variation in 
the spectral emissivity index between ~ 250 fim and ~ 5 mm. At 
these wavelengths y8 ~ 1.5 — 1.6 for both of the amC and aSil 
materials. 

The IR extinction (Fig. [5]) is dominated by the contribu- 
tion of the LamC and aSil grains. The 9.7 and 18/mi bands 
are produced by silicates while the 1 - lO/vm continuum is due 
to amorphous carbon as seen in Figs. [5] and I A. 1 1 In the DL07 
model, these two spectral components are also produced by 
amorphous silicate and big carbonaceous (graphite) grains. The 
217 nm extinction bump (see Fig.|4]i is entirel y due to aromatic 
7r — » 7r* transitions in PAHs, as proposed by llati et alJ (l2008h 
and ICecchi-Pestellini et alJ d2008[) and supported by the exper- 
imental results of ISteglich et al.l ([2010). The far-UV extinction 
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rise is also mostly produced by PAHs. This latter feature is com- 
mon to most carbonaceous grains: it is the red wing of an absorp- 
tion p eak around 750 A due to cr — > <x* transitions dJager et alJ 
2009). In our model, the bump s trength and the PAH em ission 
are thus correlated. We note that Boulange r et alJ d 19941) found 
a trend between the bump area and the IRAS 12 and 25 fim 
emission, suggesting that small amor phous carb on grains may 
also contribute to the bump (see also iMafhiil 1 9941) . The 217 nm 
and 750 A broad bands are quite weak for the BE amC sample 
that we used. We emphasize that using a carbonaceous material 
with a strong 2 17 nm band (e.g. graphite) would have resulted in 
an overestimation of this feature in our model, given the abun- 
dances required for both PAH and the small carbonaceous grains 
to reproduce the Aromatic Infrared Bands (AIB) and mid-IR 
continuum. Indeed, in the DL07 model, PAHs account for 4.6% 
of the total dust mass compared to 7.7% for our model) which 
allows for about half of the 217 nm bump to be produced by 
small graphite particles. This directly translates into a difference 
of a factor of ~ 2 in the AIB emission intensity between the two 
models as shown on Fig(6] 



The dust in our model, excited by the MMP83 ISRF, radi- 
ates a total power of 4tt j vI v dv/N H = 5.0 x 10~ 24 erg s" 1 H 1 
with a fractional contribution of 29, 7, 38 and 26% for PAH, 
SamC, LamC and aSil, respectively. This value for the total 
emitted power by dust is intermediate between those of DBP90 
and DL07 (U=l, q PAH = 4.6%), 6.9 x 10~ 24 erg s" 1 H 1 and 
4.5 x 10~ 24 erg s H , respectively, while the three models have 
about the same extinction curve (the observed one). We note that 
in the case of DBP90, the reference SED (to which the model is 
fit) was higher as can be seen in Fig|6]implying a higher dust ab- 
sorption efficiency. Indeed, we see in the bottom panel of Fig. [4] 
that the DPB90 model has the lowest albedo followed by our 
model, while the DL07 model has the highest albedo, thus ex- 
plaining the differences in emitted power noted above. As al- 
ready discussed above, the difference between our model and 
the DL07 modeled SEDs mostly resides in the AIB spectrum 
intensity while the two models agree rather well at A <: 60 fim. 



The total dust-to-gas mass ratio in our model is Mdust/Mn ~ 
10.2 x 10~ 3 , corresponding to M gas /M dus , ~ 133 (with 
Mgas/Mjj = 1 .36). This value is equal to the value for the DL07 
model (M dust /M H ~ 10.4 x 10~ 3 ) and higher than the DBP90 
model (Mdust/MH of 7.3 x 10~ 3 ). The required elemental abun- 
dances are in agreement with measurements (see TableO, ex- 
cept for iron for which our model requires ~ 30% more than 
the values inferred from both F,G stars and the Sun. The car- 
bon is shared between PAH, SamC and LamC with 32, 7 and 
61% in each respectively. Our model requires about the same 
amount of carbonaceous material and silicate as the DL07 model 
(234 ppm for C, compared to 200 ppm for our model, and 44 ppm 
of MgFeSiC»4 for DL07 compared to 45 ppm for our model). In 
our model, silicate dust represents ~ 76% (~73% in DL07) of 
the total dust mass. This mass ratio depends on the assumed den- 
sities then on the assumed form for both carbonaceous and sili- 
cate material. On the other hand, the dust material is assumed to 
be compact to keep the model as simple as possible but consid- 
ering porosity would induce a change in the carbon/silicate mass 
ratio (mainly due to change of optical properties) and decreases 
the required abundances. 
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Fig. 8. The predicted photometric band ratios Ri (mid-IR) and 
R2 (far-IR) as a function of U and the ratio between the mass 
abundance of small (Ysc) and large (Ilg) grains (see text). Note 
that the band ratios are shown on logarithmic scales. Y$g I Ylg 
is given relative to the diffuse interstellar medium value of 0.10 
for the DHGL value, which is shown by the open diamond. 



5. Tracing the nature and evolution of dust using a 
global SED 

Thanks to the Herschel and Planck missions we will soon have at 
our disposal the global dust SED, by combining these data with 
the already available near- and mid-IR data at comparable angu- 
lar resolution (IRAS, ISO and Spitzer). The availability of this 
global dust SED will allow us, for the first time, to make a sys- 
tematic study of dust evolution throughout the ISM. Indeed, dust 
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Fig. 9. The fractional contributions of the dust populations in our model to the IRAS and Planck photometric bands as a function 
of U. The bold lines emphasize the contribution of the large grains, LamC and aSil. The dust properties are those for DHGL model 
(see § l4~2b . 



evolution processes (e.g. coagulation, fragmentation) primarily 
affect the dust size distribution and the abundance ratios of the 
various dust components and are therefore directly reflected in 
the observed dust SED, which is also affected by the intensity of 
the ISRF. We first examine the effects of the latter on the model 
dust emission and the behavior of the associated fluxes in the 
photometric bands. For this purpose we scale the MMP83 ISRF 
by a factor, U, using the same nomenclature as DL07. Note that 
as we scale the entire exciting ISRF, we have U = Go, where 
Go is the ratio of the intensity of the radiation field integrated 
between 6 and 13.6 eV relative to the standard ISRF. 

Fig. [7] shows the dust SED for ISRF intensities ranging from 
U - 10~' to 10 4 . We overlay the corresponding photometric 
points for several instruments and have separated the high (left 
panel) and low (right panel) spatial resolution instruments. We 
do not show the MIPS 70 and 160yum bands that fall at the same 
wavelengths as the Herschel/PACS bands. To illustrate changes 
in the SED shape, the emission flux has been divided by U. 
As expected, we see that the emission of PAHs, which undergo 
temperature fluctuations, scales linearly with U. By contrast, the 
emission from the large grains (at thermal equilibrium with the 
ISRF) shows a more complex behavior, scaling but also shift- 
ing to higher frequencies as U increases. As U increases smaller 
and smaller grains reach thermal equilibrium with the ISRF (see 
for instance, Fig.4 of DL07) and consequently follow the same 
type of behavior as the largest grains. This is the reason why the 



SamC population shows intermediate behavior in the presented 
U range. The general behavior of our model regarding U is sim- 
ilar to that of DL07. However, a significant difference arises at 
U > 10 3 due to the fact that we use amorphous carbon rather 
than graphite. At these values of U our large grain (LamC and 
aSil) spectrum is remarkably broader, because we lack the broad 
feature in the absorption (i.e. emissivity) efficiency of graphite 
at A ~ 30 jum introduced by DL07. 

Here we illustrate how the effects of the intensity of the ISRF 
can be disentangled from effects of the relative abundances of 
small-to-large grains (i.e., PAH+SamC to LamC+aSil). The lat- 
ter variation is equivalent to a changing size distribution, and is 
to be expected when small grains coagulate onto, or disaggre- 
gate from, larger grains. We take Y, to be the mass abundances 
of dust population i (see Tab|2| and define Y$c - Ypah + Ysamc 
and Yic = YhamC + Y a su, the mass abundances of small grains 
(SG) and large grains (LG), respectively. Note, that while vary- 
ing Ysc, we keep Yuq constant and, for simplicity, assume that 
the opacity of the large grains is unchanged even though it 
should change when coagulation occurfl We also note that, even 
though the PAH/SamC abundance ratio is taken to be constant 
for this illustrative case, the relative abundance between the PAH 
band emitters and the mid-IR continuum emitters (SamC in our 



6 T his process has been shown to alter the sub -mm opacity by factors 
2 to 3 (Step nik et al.ll2003l : lParadis et all2009ah . 



10 



M. Compiegne et al.: The global dust SED: Tracing the nature and evolution of dust with DustEM 



1.0 
0.8 

0.6 

0.4 

0.2 
0.0 



SPIRE 250//m 


SPIRE 350/im 


'| Ml "Ml| 1 1 1 Mill 1 lllllll| 1 IIIIIII 

'. SPIRE 500/xm : 


^ 

PAH — — — 

- - - - SamC 

LomC 

" aSil 

LamC + aSil (LG) 




r..-^— 

- 


— i — i i-i-iiiri — nmint "T iTh hTiI - -i- 1 1 1 ntil " - r i "n iiii 

i i 1 1 1 1 1 ii i i 1 1 1 1 hi i i 1 1 1 1 hi i i 1 1 1 1 1 ii i i 1 1 1 1 1 1 


- h itinif ri m i fit" - h -Pi MTil ~ "i" i "i i f ilit "i "hum 

i i 1 1 1 1 hi i i 1 1 1 1 hi i i 1 1 1 1 1 ii i i i 1 1 1 hi i iiiiiii 


- i~i iiiiiih n-rrnfit'-TiTifi]il"i Tiiniil — 1 1 Mini 


PACS 70/im 


r PACS 10(ty*m \ 


PACS 160/xm 


- / ■ ^ 

/ / ^ 


^ -»» — 


• 

• — , 


y .x- ^ ^ 

/\ ..." 

— ■■ • ""x 


■ ..•**" 
.••"** 


• — . ■ • • 

..«- ' ^ 

• ""**' ' ■ ■ 




s 

. .1. j j mill ..ti .11 mil . .1. i J mill . .1. j .1 mill • ■!■ nirin 


n mnrf- n TTiifit' i in lull . .i. ujmil ■ -i- 1 -n mi 



10~ 1 1 10 10 2 10 3 10 4 

u 

Fig. 10. Same as Fig.|9]but for Herschel SPIRE and PACS bands. 



case) is known to vary (e .g lBerne et al.ll2007t ICompiegn e et alJ 
2008; iFlagev et alj|2009f) . Fig.g] shows the photometric ratios 
Ri a = IRAC5.S/PACS 100, R 2a = PACS 70/ S PIRE 350 (top 
panel) and R lb = IRAS 12/ IRAS 100, R 2b = IRAS 60/ HFI 350 
(lower panel) as a function of U and Y$g I Ylg- The photometric 
fluxes were obtained by integrating the spectra over the relevant 
filter transmission taking into account the color corrections. The 
ratios R\ and R 2 represent roughly the same wavelengths in case 
(a) and (b) and so consequently show a similar behavior. We see 
that for 10 3 < Y SC /Y LG < 10 and 10" 1 < U < 10 4 , the band 
ratios Ri (a and b) scale with Fsc/itc because the SG emission 
becomes dominant at small wavelengths. Conversely, the R 2 ra- 
tios (a and b) are mostly sensitive to variations in U because the 
photometric bands involved are dominated by the LG contribu- 
tion (see Fig.l9land[T0b. In this range of U and Y$gI Ylg parame- 
ter space the variations induce well discriminated differences in 
the photometric band ratios /?;. From the observed values of R\ 
and R 2 , U and Ysg/Ylg can thus be constrained for a given dust 
model. However, we note that for high U and low Ysg/Ylg Rib 
and R 2 h this discrimination is not possible because the spectrum 
is dominated by BGs down to wavelengths as short as 12 jjm. 
This degeneracy is avoided for higher U when using the shorter 
wavelength to trace SG emission, as seen with R\ a . 

In general, a photometric point measured on the dust SED 
is a sum of the contributions of the different dust components 
(composition and sizes) and the analysis of photometric data 
therefore requires prior knowledge of these contributions. In 
view of the differing behaviors of the SG and LG dust emissions, 



with respect to U (Fig.|7]i, we show in Fig.|9]and[T0]the fractional 
contributions of the model dust populations (for the DHGL prop- 
erties) in several photometric bands, as a function of U. For 
10 1 < U < 10 4 , in the IRAS 25fim band (Fig.EJ, the LG con- 
tribution gradually increases as U increases because these grains 
become warmer, as can be seen in Fig. [7] The LamC contribution 
is dominant (with respect to aSil) at A ^ 60 fxm because these 
grains are warmer, a consequence of their higher absorptivity in 
the near-IR (Fig.|ATJi. At A ~ 70 - 100 jum {A ~ 160 - 250yum), 
and for this same reason, LamC dominates at U ^ 10 2 (U ^ 5) 
while aSil does only for higher U since the LamC grains then 
emit at shorter wavelengths. LamC and aSil have about the same 
contribution at A ^ 350yum (as seen for U = 1 in Fig|2|i. While 
the wavelength of the LG emission peak moves to higher fre- 
quencies when U increases (U <: 10 2 - 10 3 ), the cold part of the 
temperature fluctuations in PAHs contributes to the emission in 
the HFI 1.4 mm (217 GHz) band (Fig.©. We note that at these 
frequencies recent theoretical work indicate that rotational emis- 
sion of PAHs for Go ^ 10 3 and gas dens i ties ~ 10 5 cm" 3 
might also contribute dHoang et al.l l2Q10h [Siisbee et alj|2010l) . 
We emphasize the fact that the dust contributions at 60 and 
100//m allow for a better determination of the LG SED around 
its emission peak and hence better mean temperature estimates 
for the LGs, a necessary step in the extraction of the thermal dust 
component in Planck data. Finally, and as discussed earlier, the 
emission of aSil grains is exp ected to be polarized u nlike that of 
the LamC grains and PAHs ( Sironi & Draine 200%. A determi- 
nation of the respective contributions of the different LGs and of 
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PAHs to the dust SED will therefore be crucial to the analysis of 
data from the Planck-HFI polarized channels (at A > 850 ymi). 

6. Summary 

The emission from large interstellar grains (a ~ 0.1 yum) is cur- 
rently being observed by the instruments onboard the Planck 
and Herschel satellites. The properties of the bulk of this dust 
are likely determined by dust evolution throughout the lifecy- 
cle of interstellar matter and are probably not just a reflection 
of the composition of unaltered 'Stardust'. The primary effect of 
dust evolution is to redistribute dust mass during growth or frag- 
mentation episodes, resulting in changes in the abundance ra- 
tio of small-to-large grains, YsgIYlg, (where small grains have 
a Z> 10 nm). Such changes are reflected in the global (IR to mm) 
dust SED, which will soon available for large parts of the sky, 
thanks to the combination of Spitzer, IRAS, Herschel and Planck 
data. Understanding dust evolution is critical for the physics and 
chemistry of the ISM but also for the optimum subtraction of the 
thermal dust contribution from the Planck data in CMB studies. 

Dust evolves in the ISM and thus shows regional variations 
in its intrinsic properties, structures and sizes. In this work we 
describe a framework dust model and provide the tools neces- 
sary for the study of dust evolution in the ISM. We present a ver- 
satile numerical tool, DustEM , that allows us to predict the dust 
SED and extinction, given the dust properties, and that is capable 
of handling a large diversity of dust types, using a knowledge of 
their size-dependent opacities and heat capacities. The design of 
DustEM allows the user to add new grain physics in a simple way. 
To enable an assessment of the impact of dust evolution we have 
defined a reference dust model for the DHGL medium, where 
the dust populations have certainly undergone growth and frag- 
mentation processing and where the calculated dust SED and ex- 
tinction represent the 'average' or 'equilibrium' observables for 
evolved dust. Our dust model comprises (i) Polycyclic Aromatic 
Hydrocarbons, (ii) amorphous carbons and (iii) amorphous sili- 
cates. We argue that hydrogenated amorphous carbon grains are 
the natural product of dust evolution in the ISM, as are amor- 
phous silicates. Although composite grains may naturally arise 
in the ISM, we did not include them here, on the basis of the re- 
sults from recent spectropolarimetric observations. The Planck 
polarized channel data will soon make available more constraints 
in this issue. 

The dust SED reflects changes in the abundances and sizes of 
grain populations and also variations in the intensity of the excit- 
ing radiation field. Using DustEM we have shown that well cho- 
sen photometric band ratios can disentangle the influence of the 
intensity of the exciting radiation field from that of changes in 
dust abundances as indicated by YsgIYlg- We show the contri- 
butions of the small and large grains in the IRAS and Planck-HFI 
channels (and similarly to Herschel channels). Such information 
allows for a reliable determination of the temperature of large 
grains, a necessary step in the quantitative subtraction of ther- 
mal dust emission from the Planck data. These results assume 
fixed dust optical properties and size distributions. However, in 
a realistic treatment of dust evolution such quantities will change 
with the local physical conditions and affect the FIR to sub-mm 
grain emissivity. Future work will require the use of IR, Herschel 
and Planck data to assess variations of Y$gI Ylg m interstellar re- 
gions where small grains may coagulate onto, or result from the 
fragmentation of, larger grains, resulting in changes in the sub- 
mm opacity of dust. Combined with the constraints from Planck 
polarized channel data, these full SED dust data should provide 
plausible constraints on the candidate materials, and their phys- 



ical properties, that can be invoked for the evolved dust in the 
ISM. 
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Appendix A: Dust properties 

Here we discuss the optical and thermal properties of the grain 
types adopted in this work. 

A.1. PAHs 

Although a global picture of the contributio ns of different PAH 
famili es to the mid-IR bands is emerging (e.g. lBauschlicher et al.l 
(2009)), it does not yet quantitatively describe how the band pro- 
files vary as a function of the physical conditions (radiation field 
intensity and gas density). As in the case of silicates (see here- 
after in the Appendix), an empirical approach can be followed 
where the r atios and profiles of the IR bands are derive d from ob- 
servations dVerstraete et alJl200U iLi & Draindl2001al) . This ap- 
proach is however degenerate, given the fact that the band ratios 
depe nd on s i ze, ion ization and hydrogenation states. Recently, 
iDraine & Lil (I2007I) (hereafter DL07) added the c onstraints pro- 
vided by the extensive theoretical database of iMalloci et al.l 
(I2007I) to define a new set of band parameters, both for neu- 
tral and cationic PAHs. Here we use the cross-sections defined 
in DL07 with the following modifications. To allow for a better 
match to the ISOCAM/ CVF spectrum of th e diffuse interstel- 
lar medium, presented in Flagey et al. (2006) and in Fig. [2] we 
increase the integrated cross-sections cr in , of the 8.6, 11.3 and 
12.7 /iin bands by 50%, 10% and 50%, respectively. We have 
not included the near-IR bands between 1 and 2 /im, which have 
litt le effect on the emission. W e added the far-IR bands defined 
by lYsard & Verstraetd (2010) from the Malloci database. The 
PAHs are assumed to be disk-like and their size, a, is defined 
to be that of an equivalent sphere of the same mass, assuming 
a density of p — 2.24 g/cm 3 . The relation between a and Nc 
is Nc — 470 (a/lnm). We assume the PAHs to be fully hydro- 
genated, with [H/C] varying as a function of Nc as defined in 
DL07. Following DL07, the same UV-visible absorption cross- 
section is assumed for neutral and charged PAHs while the same 
absorption cross-section is assumed for cations and anions. The 
size dependence of the ionization fraction of PAHs is also ob- 
tai ned in the same way as by DL0 7: it has been estimated using 
the IWeingartner & Drain e (2001b ) formalism and fo r the phys- 
ical conditions summarized in ILi & Draind (1200 lal) (see their 
Table 2) fo r the DHGL medi um. The heat capacity has been 
taken from lDraine & LH (120011) . 

A.2. Amorphous carbon 

Amorphous carbons can assume many different structures and 
compositions, depending on the way in which the sp 2 and sp 3 
carbon atoms combine wi th one another and with other species, 
especially with hydrogen. Zubko et al. (1996) compiled the op- 
tical properties of several solid hydrocarbons over a broad spec- 
tral range, from the millimeter to the EUV range. These opti- 
cal properties were d erived from laboratory measurements of 
the opacity made by IColangeli et al.l (1 19951) . In this work we 
adopt the optical properties of the BE sample: formed by ben- 
zene combustion in air, this sample is representative of rather 
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Fig.A.l. The refractive indices (m = n + ik) of amorphous carbon and amorphous silicates (left panel) and the corresponding 
absorption efficiencies obtained from a Mie calculation (right panel). 



sp 2 -rich hydrogenated amorphous carbon. The refractive index 
for the BE material is displayed in Fig. lA.ll (left). The absorp- 
tion and extinction efficiencies were obtained from Mie calcu- 
lations. Fig. lA. II (right) shows the absorption efficiencies, Q a bs, 
fo r grains of radius 1 and 100 nm. The refractive index given 
bv lZubko et alJ dl996h goes down to ~1900jt/m. We extrapolated 
the optical efficiencies down to 10 cm using the average emis- 
sivity index between 800 and 1900yum, using /3 = 1.55. This 
ft value does not depend on the siz e (for a 5= 3 ymi). For the 
heat capacitv rMichelsen et alJ (120081) argue that the graphite heat 
capacity is a good surrogate for soots and for amorphous car- 
bons: we therefore use the graphite heat capacity as defined in 
iDraine & Lil d200lh . Finally, we adopt a mass density p = 1.81 
g crrT 3 , typical of hydrogenated amorphous carbons. 

A.3. Astronomical silicate 

The observed profiles of the 9.7 and 18 /jm absorption bands in- 
dicate that the interstellar amorphous silicates are of olivine-type 
(Mg v Fe v Si04). The laboratory band profiles are, however, too 
narrow to explain the observations, suggesting that some kind 
of processing is at w ork and that other t ypes of silicates con- 
tribute to these bands (Demyk et al. 2004). Draine & Led (11984 
followed an empirical approach and defined astronomical sili- 
cates whose optical pro perties in the mid-IR were constructed 
from observational data. iLi & Draind (|2001a) further modified 
the imaginary part of the dielectric function for A > 250 fim to 
better match the high galactic latitude dust emission mea sured 
by FIR AS. Using the optical properties described in iDraind 
(2003bfl we derived the absorption and extinction efficiencies 
using Mie calculations. The right part of Fig. lA. ll shows the cor- 
responding absorption efficiencies, Qabs- Finally, we take the 



7 Available at 

http: //www. astro. princeton. edu/~draine/dust/dust. html 



heat capacity for amorphous silicates from IDraine & Lil (1200 lh 
and assume a mass density of 3.5 g crrr 3 . 

Appendix B: The temperature-dependent grain 
emissivities tool 

In this section we present a tool that allows the user to apply a 
change in the long wavelength opacity spectral index (J3) as a 
function of the wavelength and the temperature. 

The determination of dust temperatures and the dust 
spectral emissivity index, f3 (Q a bs ~ V s ), from far-IR to 
submm data is degenerate and sensitiv e to t h e errors or 
multip le dust temperatures dDupac et alJ 120031; IShettv et"al] 
2009a,b)- Nevertheless, obser vations do seem to indicate a 
tempe rature-dependence of ft dDupac et alJ 120031: iDesert et al.l 
2008), which may re flect the microscop ic properties of grains 
at low temperatures dMenv et al.l 120071) . In addition, labora- 
tory data on the sub-mm opacity of interst ellar dust analogues 
show a similar temperature dependence dAgladze et al.l [l996: 
lMennellaet"aT]|1998l; iBoudet et ai]|2005l) . Within the perspec- 
tive of the Herschel/Planck data interpretation, the possibility 
of temperature-dependent dust emissivities has been introduced 
into DustEM . This feature is implemented as a correction to the 
emissivity index (3 and takes the form: 

Qabsia, v) = Q () {a, v) ( v /v t ) S(T)H(VtM (B.l) 

where Qq is the grain emissivity without the temperature depen- 
dent index, v, is the frequency threshold below which the correc- 
tion is applied, 5(T) is the index correction and H is the thresh- 
old function. We define the correction as 5{T) = /3(T) -/3q where 
/3(T) is the temperature law for the emissivity index and /3q is 
some reference value. Such a /^-correction should then preferen- 
tially be applied to grain types whose Q a t, s behave as a power 
laws (with an index /3q) at frequencies lower than the threshold, 
otherwise fio may depend on the frequency. The /?(r)-values are 
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estimated with a user-defined, parameterized function or can be 
interpolated in a data table. The parameters, as well as the tab- 
ulated values (if used), are read from a data file. We chose the 
threshold function to be H(x) — (1 + 4stanh(x - l))/2 where s 
controls the stiffness of the transition from to 1 (Ax = x/s): for 
instance for 5 = 1 the //-values rises from to 1 for x between 
0.5 to 1.5 and the width of the transition is Ax = 1. Since the 
temperature dependence of the dust emissivity appears at long 
wavelengths (A > 30 jt/m), where the dust absorption efficiency 
is much smaller than in the visible and near-IR range (see Fig. [5), 
we apply the /^-correction only in the cooling function (emission) 
for the dust in the DustEM code. 
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